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An  integrated  mathematical  tool  is  presented  to  simultaneously  design  powertrain  and  control  strategies 
for  PEM  fuel  cell  hybrid  vehicles.  Such  an  activity  is  motivated  by  the  positive  impact  of  car  electrification 
on  transportation  related  GHG  emissions,  and  by  the  potential  of  hydrogen  as  an  intensive  energy  carrier 
to  avoid  heavy  battery  packs,  as  required  for  pure  electric  cars.  Moreover,  the  availability  of  reliable  and 
comprehensive  mathematical  tools  is  recognized  as  a  key  to  develop  cost  effective  and  highly  performing 
innovative  powertrains. 

A  database  of  hybridizing  devices’  costs  and  weights  is  developed  and  exploited  in  conjunction  with 
a  weight  model  to  size  the  powertrain  as  a  function  of  degree  of  hybridization.  In  parallel,  a  versatile 
heuristic  control  strategy  is  proposed  to  easily  adapt  control  rules  to  different  powertrains.  Then,  both 
the  weight  model  and  the  heuristic  control  strategy  are  embedded  into  an  optimization  procedure  to 
determine  the  most  convenient  PEM  fuel  cell  hybrid  vehicle  configuration. 

The  methodology  is  tested  for  a  fuel  cell  hybrid  shuttle,  whose  hydrogen  feeding  is  guaranteed  by 
a  photovoltaic  on-site  generation  system.  Such  a  case  study  not  only  is  suitable  to  assess  the  proposed 
optimization  procedure,  but  also  serves  at  indicating  the  most  promising  short-term  applications  of  fuel 
cell  hybrid  vehicles. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  the  last  years  the  increasing  interest  in  energetic  and  envi¬ 
ronmental  problems  has  given  a  strong  impulse  toward  the 
development  of  alternative  energy  conversion  devices.  Since 
transportation  systems  are  responsible  for  most  of  fossil  fuels 
consumption  and  environmental  impact,  significant  attention  has 
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been  given  to  the  development  of  alternative  powertrains.  Nowa¬ 
days,  a  large  variety  of  vehicles  are  considered  as  valid  candidates 
to  build  up  a  sustainable  mobility  worldwide.  In  the  short  term 
scenario,  hybrid  electric  vehicles  (HEV)  powered  by  internal 
combustion  engine  as  main  propulsion  device  surely  represent  the 
most  viable  solution,  thanks  to  their  mature  technology  and 
acceptable  investment  cost.  In  the  medium-term  scenario,  plugin 
HEVs,  which  are  already  ready  to  enter  the  car-market  today,  are 
expected  to  overcome  HEVs  as  best  alternative  powertrains,  espe¬ 
cially  if  the  foreseen  trends  in  Lithium-battery  technology  and 
costs,  along  with  suitable  adaptation  of  electric  grids,  will  be  met. 
The  latter  expectation  will  at  the  same  time  drive  car  makers 
toward  massive  electrification  of  their  products,  preferring  series  to 
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Nomenclature 

HT 

hydrogen  tank 

IC 

investment  cost  [€] 

A 

vehicle  frontal  area  [m2] 

ICE 

internal  combustion  engine 

AD 

annual  distance  [km] 

Jchg 

battery  charging  current 

aux 

fuel  cell  system  auxiliaries 

fdis 

battery  discharging  current 

B 

battery 

M 

mass  [kg] 

C 

cost  [€] 

m 

unit  mass  [kg  kW-1] 

c 

unit  cost  [€  kW1] 

kb 

battery  unit  cost  [€  kW-1  h-1] 

Cbus 

conventional  bus 

1<ht 

hydrogen  tank  unit  cost  [€  kW-1  h-1] 

CD 

charge  depleting 

MC 

maintenance  cost  for  PV  plant  [€] 

^Diesel 

cost  of  Diesel  [€  1_1] 

MBc 

mass  of  a  single  battery  cell  [kg] 

cE 

cost  of  the  electrolyzer  [€] 

rnHl 

hydrogen  mass  flow  [kg  s-1] 

Cr 

rolling  resistance 

NBc 

number  of  battery  cells 

cs 

charge  sustaining 

ocv 

battery  open  circuit  voltage  [V] 

cv 

conventional  vehicle 

p 

power  [kW] 

Q/vat 

cost  of  water  [€] 

p* 

nominal  power  [kW] 

Cx 

coefficient  of  drag 

PBc 

nominal  power  of  single  battery  cell  [kW] 

Eb 

battery  pack  capacity  [kW  h] 

PEM 

proton  exchange  membrane 

Em 

energy  stored  in  the  hydrogen  tank  [kW  h] 

^FCS.opt 

FCS  most  efficient  operating  point  [kW] 

EM 

electric  motor 

PHFCB 

plugin  hybrid  fuel  cell  bus 

EN 

electric  node 

Ptr 

traction  power  [kW] 

Epv 

electricity  produced  by  the  photovoltaic  plant  [kW  h] 

RB 

rule  based 

FCS 

fuel  cell  system 

^in 

battery  internal  resistance  [Ohm] 

FE 

fuel  economy  [km  l-1  or  km  kg-1] 

soc 

battery  state  of  charge 

FiT 

Feed-in  Tariff  [€  kW  h-1] 

Stack 

PEM  fuel  cell  stack 

g 

acceleration  of  gravity  [ms-2] 

T 

vehicle  transmission 

GB 

gear  box 

TAS 

total  annual  saving  [€] 

HEV 

hybrid  electric  vehicle 

fend 

driving  route  duration  [min] 

HFCV 

hybrid  fuel  cell  vehicle 

V 

vehicle  speed  [m  s-1] 

HHVh2 

hydrogen  higher  heating  value  [MJ  kg-1] 

parallel  architecture  to  maximize  battery  capacity,  with  the  final 
aim  of  extending  pure  electric  autonomy  as  much  as  possible.  Such 
a  massive  electrification  not  only  encourages  widespread  diffusion 
of  pure  electric  cars,  especially  for  urban  driving,  but  also  acts  as 
a  bridge  toward  the  introduction  of  hybrid  fuel  cell  vehicles  (HFCV), 
which  are  a  natural  evolution  of  series  plugin  HEVs  as  they  mainly 
require  to  replace  the  ICE-based  generator  by  a  much  more  efficient 
fuel  cell  system  (FCS). 

The  reason  why  HFCV  success  in  the  automotive  field  is 
considered  late  to  arrive,  despite  fuel  cell  technology  (especially 
PEM)  has  significantly  matured  in  the  last  decade,  is  strictly  related 
to  hydrogen  production,  delivery  and  on-board  storage.  One  solu¬ 
tion  that  can  reasonably  be  pursued  in  the  short  term  is  to  produce 
hydrogen  on-site  via  grid  connected  photovoltaic  panels  (PV),  thus 
making  the  production  of  hydrogen  economically  and  technically 
viable.  This  not  only  solves  the  problem  of  transporting  hydrogen, 
but  also  promotes  the  technology  and  gives  opportunities  for  early 
markets,  mainly  relying  on  currently  worldwide  available  Feed-in- 
Tariffs  (FiTs)  associated  to  PV  deployment  for  on-site  electricity 
generation  [10].  Particularly,  FiTs  can  surely  encourage  the  devel¬ 
opment  of  niche  applications  of  HFCV  powertrain,  which  of  course 
are  to  be  considered  as  highly  strategic  with  respect  to  future 
establishment  of  such  vehicular  technology. 

Moreover,  as  experienced  by  the  automotive  industry  since  the 
1980s,  both  the  appropriate  powertrain  sizing  and  definition  of 
optimal  control  strategies  have  to  be  pursued  by  developing  suited 
simulation  tools,  which  account  for  the  mutual  interferences 
between  all  the  involved  variables,  ranging  from  design  to  control 
and  exogenous  variables,  as  shown  in  Fig.  1.  Nevertheless,  while  the 
interacting  problems  of  optimal  powertrain  design  [1-3]  and  on 


board  energy  management  [4-8]  of  HFCVs  were  often  dealt  with 
separately  in  the  literature,  there  is  a  minor  contribution  on  their 
simultaneous  optimization  [9].  Therefore,  the  current  study  focuses 
on  the  development  of  an  integrated  computational  structure 
aimed  at  supporting  the  HFCV  powertrain  design,  while  guaran¬ 
teeing  that  suitable  control  strategies  are  implemented.  Such  an 
approach  is  herein  proven  suitable  to  yield  on  output  powertrain 
design  solutions  that  ensure  the  best  trade-off  among  energy 
savings,  proper  operation  of  main  hybridizing  devices  and  cost 
effectiveness  targets. 


Fig.  1.  Integrated  mathematical  tool  structure. 
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The  paper  is  structured  as  follows:  firstly,  the  modeling 
approach  followed  to  suitably  treat  the  massive  interaction  among 
design,  control,  exogenous  variables  and  constraints  is  presented; 
then,  the  versatile  heuristic  strategy  developed  to  ensure  optimal 
energy  management  be  guaranteed  on  board,  independently  of 
powertrain  configuration  and  component  specifications,  is 
described  and  discussed.  Finally,  a  case  study  is  developed  to 
highlight  the  effectiveness  of  the  proposed  methodology. 

2.  Integrated  mathematical  tool  for  optimal  design  and  on 
board  energy  management 

Fig.  1  illustrates  the  mathematical  tool  concept  that  was 
conceived  to  enhance  the  development  of  highly  performing  and 
cost-effective  HFCVs.  As  expected,  several  variable  classes  (e.g. 
design  and  control  variables)  interact  with  each  other  while 
simultaneously  defining  optimal  energy  management  and  power- 
train  sizing.  Such  interactions  are  suitably  treated  by  means  of  an 
integrated  modeling  tool,  which  was  developed  aiming  at  per¬ 
forming  optimization  analyses  in  a  wide  range  of  all  the  input 
variables  listed  in  Fig.  1.  Such  an  objective  may  only  be  pursued  via 
a  modeling  approach  that  fulfills  the  conflicting  needs  of  satisfac¬ 
tory  accuracy  and  acceptable  computational  burden,  thus  leading  to 
the  development  of  the  set  of  sub-models  described  in  the  next 
sections. 

The  whole  modeling  tool  yields  on  output  significant  metrics, 
such  as  fuel  economy,  vehicle  weight  and  cost,  which  suitably 
account  for  the  impact  of  different  design  specifications  (e.g. 
passenger  car  or  bus  application)  and  exogenous  variables  (e.g.  FiT) 
on  the  powertrain  design  optimality.  The  FiT  influence  is  of 
particular  relevance,  due  to  its  significant  role  into  the  short/ 
medium-term  deployment  of  emerging,  highly  innovative  power- 
trains,  such  as  HFCV  [11-13]. 

The  opportunity  of  evaluating  the  above  mentioned  metrics, 
taking  into  account  the  simultaneous  effect  of  different  variables 
variation,  will  in  turn  make  it  possible  to  single  out  powertrain 
designs  that  meet  the  conflicting  needs  of  cost  effectiveness  and 
energy  performance,  while  ensuring  optimal  on  board  manage¬ 
ment  of  energy  flows  be  performed.  This  latter  objective  was  ach¬ 
ieved  in  this  paper  by  developing  a  versatile,  rule-based  control 
strategy,  conceived  in  such  a  way  as  to  ensure  that  specifications- 
independent  control  strategies  are  selected  prior  to  be  fed  as 
input  to  the  integrated  tool,  as  shown  in  the  upper  left  corner  of 
Fig.  1. 


2.1.  Mass  model 

A  parametric  model  was  set  up  to  assess  the  impact  of  hybrid¬ 
ization  on  vehicle  mass.  By  referring  to  the  vehicle  architecture 
shown  in  Fig.  2,  the  mass  of  a  HFCV  can  be  obtained  by  adding  the 
mass  of  major  hybridizing  devices  to  the  vehicle  body  mass.  This 
latter  is  derived  from  the  mass  of  the  reference  conventional 
vehicle  (CV)  by  subtracting  the  contributions  due  to  the  original 
gear  box  and  ICE  from  the  original  CV  mass,  as  follows: 

Mbody  =  Mcv  -  Pice, cv  '  (micE  +  ^gb)  ( 1 ) 

Then,  vehicle  mass  can  be  determined  adding  the  hybridization 
devices  and  by  imposing  that  the  HFCV  power  to  weight  ratio  (i.e. 
pptw)  equals  CV  one,  thus  ensuring  HFCV  guarantees  the  same  CV 
acceleration  performance: 


Mhfcv  =  Mbody  +  PFCSmFCs  +  PFMmEM  +  MbcJVbc  +  MHt  (2) 


Fig.  2.  Series  architecture. 


NBc 


(3) 


MHFcv 


^EM 

Pptw 


P* 

1  EM 

^ICE,CV 


Mcv 


(4) 


Table  1  lists  the  unit  mass  here  assumed  for  each  powertrain 
component. 

It’s  worth  remarking  here  that  the  nominal  power  of  the  bat¬ 
tery’s  single  cell  (Pgc)  was  assumed  hereinafter  constantly  equal  to 
1.25  kW  [28];  therefore,  the  single  cell  mass  can  be  obtained  as 
a  function  of  either  specific  energy  or  specific  power,  which  are 
univocally  linked  to  each  other  as  shown  in  Fig.  6.  This  is  a  basic 
point  in  the  automatic  hybrid  powertrain  sizing  that  is  performed 
by  the  proposed  mathematical  tool. 


2.2.  Cost  model 

This  section  focuses  on  the  development  of  a  HFCV  cost  model, 
which  the  economical  feasibility  analysis  discussed  in  Section  4  is 


PEM  FC  efficiency 


Fig.  3.  FCS  efficiency  curves  estimated  from  measurements  performed  at  eProLab- 
UNISA. 
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Fig.  4.  Kirchoff  law-based  approach  to  battery  modeling. 


Fig.  5.  (a)  Single  Lithium  cell  OCV  vs.  SOC;  (b)  single  Lithium  cell  internal  resistances  vs.  SOC. 


based  on.  Such  a  model  allows  estimating  the  HFCV  cost,  which 
results  to  be  greater  than  conventional  vehicle  by  the  amount 
required  to  hybridize  the  powertrain.  First  of  all,  the  vehicle  body 
cost  is  determined  subtracting  the  cost  of  the  internal  combustion 
engine,  no  longer  present  in  the  FIFCV,  to  CV  global  cost,  as  follows: 

Q)ody  =  Qv  -  QCE'^ICE  (5) 

where  the  ICE  unit  cost  was  assumed  equal  to  24  €  kW-1.  There¬ 
fore,  the  FIFCV  overall  cost  can  be  estimated  taking  into  account  the 
additional  costs  due  to  hybridizing  devices’  addition: 

CHFCV  =  Q>ody  +  CFCS  ’  PFCS  +  *  EB  +  CEM  ’  PEM  +  EHT  ’  EHT  (6) 

where  the  EM  unit  cost  (including  the  inverter)  equals  16.8  €  kW-1 
[20,21  ].  On  the  other  hand,  the  unit  cost  of  fuel  cells,  hydrogen  tank 
and  batteries  is  closely  linked  to  production  volumes.  For  this 
reason,  three  cost  scenarios  were  assumed  (see  Table  2)  for  the 
main  hybridizing  devices,  aiming  at  correctly  estimating  future 
trends  in  FIFCV  economical  feasibility.  Particularly,  the  dependence 
of  PEM  FCS  cost  on  production  volume  was  derived  from  Ref.  [22], 
while  the  battery  cost  trend  proposed  in  Table  2  was  adapted  from 
the  Lithium-battery  price  forecast  proposed  in  Ref.  [23].  Instead, 
the  volume  1  cost  of  the  hydrogen  tank  is  assumed  equal  to  the 
average  cost  reported  in  Ref.  [24],  whereas  the  costs  associated  to 


volume  2  and  3  were  estimated  assuming  a  proportional  relation¬ 
ship  with  FCS  cost  trend.  It  is  worth  mentioning  here  that  the 
hydrogen  tank  cost  assumed  for  volume  3  well  agrees  with  the 
value  given  in  Ref.  [25]. 


2.3.  Longitudinal  vehicle  model 

Fuel  economies  in  this  paper  are  evaluated  by  means  of  a  back¬ 
ward  longitudinal  vehicle  model  developed  in  Matlab®  environ¬ 
ment,  whose  basic  equations  are  presented  below. 

Road  load  is  estimated  as: 

Ptr  =  MHFcvg'i'-[CrCOs(a)  +  sin(a)]  +  0.5pCxAv3  +  Mef{^v 

(7) 


Specific  power  vs.  Specific  Energy 


Table  1 

Components  unit  mass  [12—19]. 


Internal  combustion  engine  [kg  1<W_1]  2 

Gear  box  [kg  kW-1]  0.478 

Electric  motor  (including  inverter)  [kg  1<W-1]  1 

PEM  Fuel  cell  system  [kg  1<W_1]  3.70 

Hydrogen  tank  [kW  h  kg-1]  1.9 


Fig.  6.  Relationship  between  Li-ion  battery  specific  power  and  specific  energy. 
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Table  2 

Cost  trends  assumed  for  the  main  HFCV  hybridizing  devices. 


Annual  production  volumes 

Volume  1 
(1000  units) 

Volume  2 
(100,000  units) 

Volume  3 
(500,000  units) 

Fuel  cell  system  [€  kW-1] 

175 

54 

39 

Battery  [€  kW  h-1] 

1500 

1000 

500 

Hydrogen  tank  [€  kW  h-1] 

81 

25 

18 

where  a  and  p  are  the  road  grade  and  air  density,  respectively,  while 
Meff  equals  1.1  MHfcv  to  suitably  account  for  rotational  inertia.  For 
non  negative  Ptr  values,  the  mechanical  power  to  be  supplied  by  the 
EM  is  (see  Fig.  2): 

Pem  =  ^  if  Ptr  >  o  (8) 

Vt 

Pem  can  also  be  expressed  as  a  function  of  FCS  and  battery 
power,  as  follows: 

^EM  =  4em(^FCS4fCS +^b)  if  ^tr  >  0  (9) 


SOC  P  [kW] 


Time  [min] 


Fig.  7.  Qualitative  description  of  thermostatic  management  for  a  HFCV. 


where  the  p  variables  correspond  to  efficiency  terms.  On  the  other 
hand,  when  Ptr  <  0,  the  regenerative  braking  mode  is  active, 
resulting  in  the  following  expression  for  the  electrical  energy 
delivered  by  the  EM: 

Pem  =  Ptr??T??EM  if  Ptr  <  o  (10) 


good  qualitative  agreement  between  the  FCS  efficiency  curve 
acquired  at  UNISA  PEM2  test  bench  and  previously  published 
curves  [9,27].  The  values  estimated  via  Eq.  (13)  are  hereinafter  used 
to  evaluate  the  amount  of  H2  consumed  by  HFCVs  by  means  of 
Eq.  (12). 


During  regenerative  braking,  battery  can  be  charged  by  FCS  also, 
thus  the  following  equation  holds  for  negative  Ptr  values: 


2.5.  Battery  model 


Pb  =  Pem  ~  PecsVecs  if  Ptr  <  0  (11) 

It  is  worth  remarking  here  that  EM  efficiency  is  computed  by 
means  of  normalized  maps  derived  from  the  model  library 
proposed  in  Ref.  [26],  whereas  FCS  and  battery  efficiencies  are 
estimated  via  the  respective  sub-models  presented  in  Sections  2.4 
and  2.5. 

2.4.  Modeling  of  fuel  cell  system  efficiency 

In  this  paper,  a  proton  exchange  membrane  fuel  cell  was 
considered.  In  order  to  estimate  hydrogen  fuel  economy  depen¬ 
dency  on  main  HFCV  design  and  control  variables,  a  normalized 
efficiency  map  was  inferred  by  curve  fitting  the  experimental  data 
acquired  at  the  PEM  FC  test  bench  available  at  the  Energy  and 
Propulsion  Laboratory  -  University  of  Salerno  (eProLab  -  UNISA). 
Particularly,  both  gross  FC  stack  power  and  auxiliary  power 
adsorption  were  measured  to  calculate  net  FCS  efficiency,  as 
follows: 


4  FCS  = 


Pstack  ~  Paux 

mH2HHVH2 


(12) 


Then,  the  experimental  FCS  efficiency  points  were  curve  fitted 
by  the  following  rational  polynomial  function: 


The  battery  pack  model  estimates  battery  state  of  charge  (SOC) 
and  current  as  a  function  of  the  actual  electric  power  (i.e.  positive 
for  discharge  and  negative  for  charge).  The  actual  current  is 
computed  starting  from  the  electrical  power,  by  applying  the 
Kirchoffs  law  to  the  equivalent  circuit  shown  in  Fig.  4.  The 
dependence  of  main  equivalent  circuit  parameters  (i.e.  OCV  and  R[n) 
on  SOC  for  Lithium  battery  technology  was  modeled  adapting  the 
experimental  data  provided  in  Refs.  [3]  and  [28].  Such  nonlinear 
dependencies  are  illustrated  in  Fig.  5(a)  and  (b),  respectively. 
Particularly,  the  latter  figure  provides  the  R[n  values  to  be  exploited 
in  the  Kirchoff  law-based  calculation  of  single  cell  voltage,  as  shown 
in  Fig.  4. 

In  order  to  account  for  the  impact  of  battery  features  (i.e.  high 
power  or  high  energy  type)  on  HFCV  powertrain  performance  and 
cost  effectiveness,  a  further  literature  analysis  [28,29]  was  carried 
out  to  infer  the  relationship  between  battery  specific  power  and 
specific  energy,  resulting  in  the  inverse  correlation  shown  in  Fig.  6. 
Such  an  inverse  correlation,  when  coupled  with  both  nominal  FCS 
power  and  nominal  battery’s  single  cell  power,  allows  precise 
estimation  of  the  whole  battery  pack  mass  as  a  function  of  either 
specific  energy  or  specific  power. 

3.  Rule-based  energy  management  of  series  hybrid  fuel  cell 
vehicles 


-  =  Pl^FCS  +P2 

FCS  4i^fcs  +  foPpcs  + 


(13) 


Fig.  3  illustrates  the  comparison  between  experimental  and 
modeled  efficiency  values,1  highlighting  the  accuracy  achieved  by 
the  model  expressed  by  Eq.  (11).  It  is  also  worth  pointing  out  the 


The  availability  on  HFCVs  of  both  an  electric  generator  (i.e.  the 
fuel  cell  system)  and  a  battery  pack  makes  it  possible  to  operate 
HFCVs  in  plugin  mode,  as  shown  in  Fig.  7.  At  the  beginning  of  the 
driving  route,  charge  depleting  (CD)  mode  is  active,  thus  the  battery 
capacity  is  firstly  exploited  to  run  the  HFCV  powertrain  as  a  pure 
electric  vehicle  [30].  Afterward,  once  the  minimum  state  of  charge 
SOCf  has  been  reached,  the  plugin  control  strategy  switches  to  the 


1  Both  experimental  and  model  values  were  omitted  in  Fig.  3  for  confidentiality  _ 

reasons.  2  www.eprolab.unisa.it. 
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charge-sustaining  (CS)  task  (see  Fig.  7),  through  which  the  targeted 
final  state  of  charge  SOCf  is  reached  at  the  end  of  the  cycle.  Since 
HFCV  configuration  necessarily  corresponds  to  the  series  architec¬ 
ture  (see  Fig.  1 ),  the  use  of  a  simple  thermostatic  strategy  is  suitable 
to  accomplish  the  CS  control  objective  [31  ].  Thermostatic  CS  control 
of  series  hybrid  vehicles  was  actually  proven  to  be  highly  effective 
as  compared  to  ideal  fuel  economies  achievable  using  a-priori 
based  dynamic  programming  optimization  techniques  [32]. 

Fig.  7  indicates  that  the  thermostatic  control  strategy  is  adopted 
as  soon  as  the  CD  phase  expires,  i.e.  when  battery  state  of  charge 
decreases  from  fully  charge  conditions  (SOC  =  1 )  down  to  SOCf.  The 
battery  is  initially  further  depleted  until  SOC  becomes  lower  than 
a  given  threshold,  set  to  SOQ0  =  SOCf  -  dSOC.  Afterward  FCS  is 
turned  on  at  the  assigned  power  level  and  switches  off  when  the 
upper  threshold  SOCup  =  SOCf  +  dSOC  is  reached.  The  procedure  is 
repeated  until  the  end  of  the  driving  cycle,  thus  ensuring  charge- 
sustaining  operation  of  the  battery  pack  over  the  entire  CS  phase 
be  achieved.  The  simple  thermostatic  strategy  shown  in  Fig.  7  can 
be  further  improved  by  optimizing  the  parameters  Pfcs  and  dSOC  in 
such  a  way  as  to  well  adapt  on  board  energy  management  to 
different  driving  conditions  and  powertrain  configurations.  The 
following  section  describes  in  detail  the  methodology  developed  to 
optimize  such  parameters. 


3.1.  Rule-based  thermostatic  management  ofHFCVs  in  charge 
sustaining  operation 

The  optimal  values  of  Pfcs  and  dSOC  are  estimated  online  as 
a  function  of  current  traction  power  demand,  as  follows: 

Pfcs  =  /(Ptr)  (14) 

dSOC  =  g(Ptr)  (15) 

where  the  functions /and  g  correspond  to  look-up  tables  developed 
offline  extending  the  procedure  proposed  in  Refs  [32]  and  [33]  to 
series  FIFCV.  Such  a  procedure  aims  at  determining  optimal  /and  g 
functions  (i.e.  Eqs.  (14)  and  (15))  by  solving  the  following  minimi¬ 
zation  problem  over  an  extended  range  of  average  traction  power 
demand  (Ptr): 


minPFCS,AtFc-0FF  J  ™h2(Pfcs>  Atpcs— OFF)dt  (16) 

subject  to  the  constraint: 

SOC(end)  =  SOCf  (17) 

Solving  the  above  minimization  problem  corresponds  to  iden¬ 
tifying  the  best  engine  intermittency  (i.e.  AtFCs-oFF,  which  is  the 
time  during  which  FCS  is  kept  off)  and  power  level  at  which  the  FCS 
should  be  operated  for  a  given  traction  power  demand  over 
a  specified  time-horizon,  here  set  to  th  =  10  min  (see  Fig.  8). 

While  the  function  /  (i.e.  Eq.  (14))  is  directly  determined  by 
solving  the  above  described  minimization  problem  (i.e.  Eqs.  (16) 
and  (17)),  the  g  relationship  among  dSOC  and  average  traction 
power  demand  is  determined  as  shown  in  Fig.  8,  where  the  optimal 
SOC  trajectory  is  depicted:  SOC  initially  decreases  as  the  FCS  is 
imposed  to  be  initially  off;  then,  due  to  FCS  switching-on,  SOC  trend 
is  inverted  allowing  to  satisfy  the  constraint  expressed  by  Eq.  (17). 
The  halved  difference  between  SOCf  and  SOCmin  (see  Fig.  8)  is  then 
used  to  evaluate  the  SOC  excursion  (i.e.  dSOC)  to  be  adopted  in  the 
rule-based  control  strategy.  Of  course,  effective  on-line  imple¬ 
mentation  of  both  rules  (i.e.  Eqs.  (14)  and  (15))  entails  estimating 
average  traction  power  demand.  Particularly,  the  latter  variable  can 


Battery  SOC  [/] 


Fig.  8.  dSOC  indirect  determination  from  the  optimal  SOC  trajectory  addressed  by  the 
solved  minimization  problem  (see  Eqs.  (16)  and  (17))  for  an  assigned  Ptr  value. 


be  updated  every  th  time  horizon  by  means  of  either  an  a-priori  (i.e. 
Eq.  (18))  or  an  a-posteriori  (i.e.  Eq.  (19))  method: 

£;+th 

Ptr(0|,<t<,+th  =  l  J  P*m  (18) 


^  /  Mt)dt  (19) 

where  i  is  an  integer  varying  in  the  range  [0  -r-  (tend/fh)l-  Of  course 
the  a-priori  strategy  becomes  suitable  for  on  board  application 
once  coupled  either  to  GPS  information  or  model-based  forecasting 
of  average  traction  power  demand. 

Figs.  9  and  10  explain  the  variation  of  each  heuristic  rule  (see 
Eqs.  (14)  and  (15)).  The  analysis  of  both  figures  confirms  the 
conclusions  previously  drawn  by  the  authors  on  the  intrinsic 
features  of  the  proposed  rule-based  (RB)  control  strategy  for  series 
hybrid  architectures  [33].  In  the  current  paper,  the  algorithm  was 
further  refined  by  normalizing  both  rules  (i.e.  Eqs.  (14)  and  (15)), 
thus  making  them  sufficiently  versatile  to  be  successfully  deployed 


Fig.  9.  Variation  of  optimal  PFcs  (i.e.  RB  curve)  as  function  of  average  traction  power  - 
(i.e.  Eq.  (14)). 
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Fig.  10.  Variation  of  optimal  dSOC  (i.e.  RB  curve)  as  function  of  average  traction  power 
-  (i.e.  Eq.  (15)). 

over  a  wide  range  of  HFCV  designs  and  applications.  Particularly, 
both  in  Figs.  9  and  10,  the  abscissa  is  normalized.  The  effective  input 
values  of  the  /and  g  look-up  tables  (i.e.  Eqs.  (14)  and  (15))  are  ob¬ 
tained  multiplying  such  abscissa  by  PpCS.  In  such  a  way,  it  is  possible 
to  account  for  the  fact  that  the  higher  the  value  of  PpCS,  the  higher 
will  be  the  value  of  optimal  Pfcs,  as  shown  in  Fig.  3. 

Regarding  the  ordinates  of  Figs.  9  and  10,  the  de-normalization 
of  the  rule  expressed  by  Eq.  (14)  is  straightforward,  in  that  its 
output  only  needs  to  be  multiplied  by  PpCS  to  estimate  the  best 
power  level  at  which  the  FCS  should  be  operated.  On  the  other 
hand,  it  is  worth  providing  further  details  on  how  the  rule 
expressed  by  Eq.  (15)  is  obtained:  for  a  given  battery  pack  design, 
the  ordinates  of  Fig.  10,  expressed  in  kW  h,  are  divided  by  the 
overall  battery  capacity,  also  expressed  in  kW  h.  This  way,  the  rule 
illustrated  in  Fig.  10  yields  higher  dSOC  values  for  smaller  battery 
packs  and  vice  versa,  as  expected. 

3.2.  Verification  ofRB  strategy  energy  performance 

As  mentioned  above,  the  reliability  of  thermostatic  energy 
management  for  series  hybrid  vehicles  was  already  tested  against 
both  dynamic  programming  and  genetic  algorithm  off-line  opti¬ 
mization  [32,33].  In  this  section,  a  comparative  analysis  is  per¬ 
formed  to  highlight  the  effectiveness  of  the  RB  strategy  described  in 
the  previous  section  to  further  increase  energy  benefits  associated 
to  fuel  cell  hybridization  of  vehicles  as  compared  to  simple 


New  European  Driving  Cycle  (NEDC) 


Table  3 

Conventional  passenger  car  specifications. 


Length  [mm] 

4450 

Width  [mm] 

1770 

Height  [mm] 

1550 

ICE  power  [kW] 

75 

Mass  [kg] 

1300 

Coefficient  of  drag  [/] 

0.33 

Coefficient  of  rolling  friction  [/] 

0.01 

Rolling  radius  [mm] 

500 

thermostatic  management.  The  longitudinal  vehicle  model 
described  in  Section  2.3  was  applied  to  a  typical  passenger  car  (see 
specifications  in  Table  3),  which  was  considered  for  hybridization 
with  fuel  cells.  Several  HFCV  configurations  were  designed 
applying  the  models  developed  in  Section  2,  assuming  that:  i)  PpCS 
varies  in  the  range  20-60  kW;  ii)  battery  specific  energy  is 
constantly  equal  to  215  W  h  kg-1;  iii)  the  hydrogen  tank  mass  and 
energy  (see  Eqs.  (2)  and  (6))  guarantee  up  to  400  km  autonomy. 
Fig.  12  shows  the  hydrogen  fuel  economies  simulated  for  the  above 
set  of  HFCVs  on  the  new  European  driving  cycle  (see  Fig.  11).  The 
comparative  analysis  highlights  how  the  RB  energy  management 
outperforms  simple  thermostatic  management,  whose  corre¬ 
sponding  values  were  simulated  assuming  Pfcs  and  dSOC  be 
constantly  equal  to  Pfcs, opt  and  0.01,  respectively.  The  benefits  are 
particularly  evident  in  case  a-priori  estimation  of  RB  strategy  input 
(i.e.  average  traction  power)  is  performed,  thus  indicating  the  need 
for  an  on  board  predictor  of  future  average  power  required  at 
vehicle  wheels.  The  results  illustrated  in  Fig.  12  well  agree  with 
similar  analyses  available  on  the  public  domain  [9],  thus  demon¬ 
strating  the  suitability  of  the  proposed  RB  strategy  to  single  out  the 
most  energy  saving  HFCV  design  (i.e.  corresponding  to 
Ppcs  =  40  kW)  once  embedded  into  an  integrated  mathematical 
tool  as  the  one  proposed  in  this  paper.  It  is  worth  remarking  here 
that  the  proposed  RB  management  strategy  does  not  need  to  be 
optimized  as  a  function  of  powertrain  degree  of  hybridization,  thus 
highlighting  once  again  its  contributing  value  with  respect  to  the 
methodology  previously  developed  by  Kim  and  Peng  [9]. 

Fig.  13  shows  the  Pfcs  trajectory  simulated  for  the  most  energy 
saving  solution,  both  when  enabling  RB  and  simple  thermostatic 
strategy.  As  expected,  and  as  more  clearly  highlighted  in  the  close 
window  shown  in  Fig.  13(b),  the  RB  strategy  allows  avoiding  battery 
overcharging  during  low  power  demand  phases,  thus  resulting  in 


( - 1  PFCS=PFCS,opt  RB  a-priori 

RB  a-posteriori  ---Poli.  (RB  a-priori) 


FCS  nominal  power  [kW] 

Fig.  12.  Comparison  of  fuel  economies  simulated  with  RB  and  simple  thermostatic 
energy  management. 
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NEDC  power  demand 


Time  [s] 


Fig.  13.  (a)  Comparison  of  simulated  fuel  cell  power  trajectories  for  the  best  fuel  economy  powertrain  (i.e.  corresponding  to  PpCS  =  40  kW,  as  shown  in  Fig.  12).  (b)  Close  window 
of  (a). 


Table  4 

Characteristics  of  the  conventional  Diesel  bus. 


Length  [mm] 

9420 

Width  [mm] 

2330 

Height  [mm] 

2885 

Power  [kW] 

160 

Mass  [kg] 

14,655 

Coefficient  of  drag  [/] 

0.5 

Coefficient  of  rolling  friction  [/] 

0.012 

Rolling  radius  [mm] 

500 

Diesel  fuel  economy  [km  l-1] 

2.1 

a  more  efficient  use  of  the  battery  pack.  Moreover,  Fig.  13(a)  indi¬ 
cates  how  the  RB  strategy  forces  the  FCS  to  provide  boost  power 
during  peak  power  demand,  thus  maximizing  direct  flow  of  energy 
from  FCS  to  the  wheels,  which  in  turn  reduces  charging/discharging 
dissipations  associated  to  excessive  energy  storage  into  the  battery 
pack. 


4.  Case  study 

The  case  study  focuses  on  a  real  world  potential  early  market  for 
HFCV,  namely  the  development  of  a  plugin  hybrid  fuel  cell  bus 
(PFIFCB)  to  be  used  in  place  of  the  conventional  bus  (Cbus,  see 
specifications  in  Table  4)  that  is  currently  used  as  a  shuttle  at  UNISA 
campus.  Hydrogen  feeding  and  battery  recharging  are  guaranteed 
by  an  on-site  PV  generation  system.  Such  a  solution  allows,  on  one 
hand,  to  deeply  exploit  FiT  associated  to  photovoltaic  energy 
generation  and,  on  the  other,  to  produce  clean  energy  and  thus 
develop  a  vehicle  with  no  environmental  impact.  Table  5  summa¬ 
rizes  the  main  features  of  the  case  study.  Furthermore,  the  latter 
table  specifies  the  value  of  several  variables  introduced  in  the 
following  Eqs.  (20)— (22). 

In  order  to  assess  the  economic  feasibility  of  replacing  the 
existing  bus  by  a  fuel  cell  vehicle,  a  pay-back  analysis  was  per¬ 
formed.  In  this  analysis  the  total  investment  is  considered  equal  to 

Table  5 

Characteristics  of  the  case  study. 


Daily  distance  [km  day-1]  156.7 

Working  days  [/]  238 

Annual  distance  [km  year-1]  37,295 

Cost  of  Diesel  fuel  [€  l1]  1.80 

Cost  of  water  [€  m-3]  0.75 

Cost  of  PV  panels  [€  1<W-1]  3270 

FiTa  [€  kWh-1]  0.233 

PV  energy  per  peak  power  [kW  h  1<W-1]  1300 


a  The  assumed  value  corresponds  to  2011  Italian  FiT  for  PV 
diffusion  enhancement  [34]. 


the  difference  between  the  new  PHFCB  and  the  current  conven¬ 
tional  bus  cost,  which  is  then  added  to  the  photovoltaic  plant  and 
electrolyzer  cost: 


IC  =  CpHFCB  -  Qlbus  +  Cpv'Ppv  +  (20) 

It  is  worth  pointing  out  here  that  both  unit  cost  and  nominal 
power  of  photovoltaic  plant  [35]  and  electrolyzer  [36]  depend  on  the 
PHFCB  characteristics  and,  in  turn,  on  its  energetic  performance  (i.e. 
hydrogen  fuel  economy  FEH2  [km  kg-1]),  to  be  simulated  as  a  func¬ 
tion  of  powertrain  specifications  by  means  of  the  flexible  mathe¬ 
matical  tool  presented  and  discussed  in  the  previous  sections. 
Particularly,  the  nominal  PV  power  to  be  installed  to  guarantee  pure 
photovoltaic  hydrogen  generation  is  computed  as  follows: 


P;V[1<W] 


AD  HHVh.-KH-SSOO-1  \  1 

f%  %  +  py 


(21) 


where  hE  and  EPV  Ppv-1  are  the  electrolyzer  efficiency  and  the  PV 
energy  [kW  h]  per  peak  power  [kW],  respectively,  while  Epi  is  the 
annual  electrical  energy  consumed  to  compensate  for  plugin 
operation  of  the  vehicle  (i.e.  to  recharge  the  battery  at  the  end  of 
the  day  to  increase  SOC  back  to  fully  charge  condition).  It  is  worth 
remarking  here  that  the  value  assumed  for  t]e  in  this  work  (i.e.  0.63 
and  0.53  if  higher  or  lower  H2  heating  value  is  assumed,  respec¬ 
tively  [36])  is  safe  enough  to  be  representative  of  the  entire 
hydrogen  production  chain  via  alkaline  electrolysis,  including  the 
energy  spent  for  pressurization  [37].  Eq.  (21)  also  highlights  that 
the  only  cost  associated  to  the  electricity  required  for  hydrogen 
production  corresponds  to  PV  investment  cost. 
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Fig.  15.  Parametric  analysis  outcomes:  (a)  volume  1;  (b)  volume  2;  (c)  volume  3. 


Table  6 

Simulated  performance  for  the  best  case  shown  in  Fig.  15  (i.e.  nominal  fuel  cell 
system  power  =  30  kW,  battery  specific  energy=110  kW  h). 


Mass  [kg] 

15,571 

FCS  nominal  power  [kW] 

30 

Battery  pack  capacity  [kW  h] 

28.4 

Number  of  battery  cells  [/] 

106 

EM  nominal  power  [kW] 

161.5 

Hydrogen  fuel  economy  in  CS  phase  [km  kg-1] 

20.72 

Diesel  equivalent  fuel  economy  [km  l-1] 

4.21 

Energy  consumption  in  CD  phase  [W  h  km-1] 

745 

Regarding  electrolyzer  cost,  in  this  work  it  is  evaluated  as 
a  function  of  daily  hydrogen  production,  by  means  of  the  rela¬ 
tionship  proposed  in  Ref.  [36].  It  is  also  worth  mentioning  that  such 
relationship,  which  holds  valid  as  it  is  for  the  volume  1  scenario 
(see  Table  2),  is  adapted  to  production  volumes  2  and  3  by  pro¬ 
portionally  scaling  its  output  with  respect  to  FCS  cost  reduction. 

The  total  annual  savings  is  instead  related  to  the  saved  Diesel 
consumption  and  FiT,  to  which  however  the  cost  of  water  needed 
for  electrolysis  and  the  PV  annual  maintenance  cost  (here  set  equal 
to  1%  of  entire  PV  plant  cost)  must  be  subtracted: 

TAS  =  (=£5_)  .  cDiesel  +  EPV  •  FiT  -  Cwat  -  MC  (22) 

VhtCbus/ 

The  payback  can  then  be  evaluated  as: 

Payback  =  ^  (23) 

The  pay-back  analysis  was  applied  to  the  three  production 
volumes  shown  in  Table  2,  thus  accounting  for  the  impact  of  market 
penetration  of  hybridizing  devices  on  PFIFCB  comprehensive  cost. 
In  addition,  for  each  production  volume  a  parametric  analysis  was 
performed  to  evaluate  payback  variation  as  a  function  of  nominal 
fuel  cell  system  power  and  battery  specific  energy,  in  the  ranges 


[30, 130]  kW  and  [110, 220]  kW  h,  respectively.  Such  an  analysis  was 
conducted  simulating  the  PFIFCB  on  a  real  world  bus  driving 
schedule,  namely  the  London  bus  cycle  (see  Fig.  14).  In  this  case,  the 
mass  and  cost  models  (see  Sections  2.1  And  2.2)  were  deployed 
assuming  that  hydrogen  tank  size  ensures  up  to  150  km  autonomy, 
which  is  consistent  with  the  case  study  characteristics  (see  Table  5) 
and  with  the  fact  that  PFIFCB  is  run  as  a  plugin  vehicle.  Battery  mass 
was  computed  by  suitably  coupling  Eqs.  (2)— (4)  to  the  information 
on  battery  specific  power  retrievable  from  Fig.  6.  It  is  also  worth 
remarking  here  that  the  versatile  features  of  the  RB  management 
strategy  presented  in  Section  3.1  were  exploited  by  adapting  the 
rules  described  in  Figs.  9  and  10  to  the  varying  values  of  FCS 
nominal  power  and  battery  specific  energy  considered  in  the 
parametric  analysis. 

Fig.  15  shows  the  results  of  the  above  described  analysis.  As 
expected,  projected  reduction  in  hybridizing  devices  cost  will 
enhance  economic  feasibility  of  the  proposed  solution.  More 
interesting  to  note  is  the  fact  that  a  range-extender  like  powertrain 
design  (i.e.  low  PpCS)  emerges  as  the  best  compromise  between 
energetic  performance  and  cost  effectiveness  (see  Table  6).  There 
are  two  main  justifications  of  such  a  result.  First  of  all,  buses  do  not 
suffer  so  significantly  from  hybridization  mass  increase,  thus 
resulting  in  negligible  differences  in  energetic  performance  when 
switching  from  one  degree  of  hybridization  to  another  (i.e.  by  either 
reducing  or  increasing  nominal  FCS  power).  Moreover,  reducing 
nominal  FCS  power  results  in  higher  battery  nominal  power  (see 
Section  2.1),  which  in  turn  causes  battery  storage  capacity  to 
increase  without  needing  high  energy  density  (see  Fig.  6).  This  way, 
plugin  effect  is  better  exploited  as  compared  to  high  PpCS  solutions 
(see  Fig.  15),  without  requiring  such  a  high  investment  cost  for  both 
battery  pack  and  fuel  cell  system.  In  this  context,  the  availability  of 
a  previously  optimized  control  strategy  (see  Section  3),  which  is 
suitable  to  be  embedded  within  the  integrated  mathematical  tool 
proposed  in  this  paper,  was  shown  to  be  effective  in  ensuring 
optimization  analyses  be  conducted  without  over-simplifying 


Time  [min]  Time  [min] 


Fig.  16.  (a)  Comparison  of  simulated  SOC  trajectories,  (b)  Close  window  of  (a). 


M.  Sorrentino  et  al.  /  Journal  of  Power  Sources  221  (2013)  308-317 


317 


hypotheses.  Usually,  one  common  over-simplification  relates  to 
how  is  accounted  for  the  dependence  of  powertrain  energetic 
performance  on  design  specifications  and  vehicle  application  (i.e. 
prevailing  urban  use  of  the  bus-shuttle,  which  the  case  study 
focuses  on). 

Finally,  the  selection  of  a  range-extender  like  PHFCB  design  appears 
satisfactory  even  from  a  battery  lifetime  point  of  view,  as  it  emerges 
from  Fig.  16,  where  the  SOC  trajectories  simulated  for  both  the  best  and 
the  PpCS  =  110  kW  solutions  are  compared.  Particularly,  the  SOC 
trajectory  associated  to  PpCS  =  110  kW  exhibits  a  higher  frequency 
oscillation  around  the  targeted  SOCf,  which  in  turn  might  lead  to 
higher  number  of  equivalent  DOD  with  respect  to  PpCS  =  30  kW. 

5.  Conclusions 

An  integrated  mathematical  tool  was  conceived  to  simulta¬ 
neously  optimize  powertrain  sizing  in  hybrid  fuel  cell  vehicles  and 
their  control  strategies.  The  models  were  developed  aiming  at 
finding  the  best  compromise  between  conflicting  needs,  such  as 
accuracy  and  computational  burden,  thus  making  the  proposed 
methodology  suitable  to  run  extended  optimal  design  oriented 
analyses. 

Moreover,  the  integration  of  the  modeling  tool  with  a  versatile 
heuristic  strategy  was  proven  to  be  effective  to  guarantee  high 
performance  be  achieved  independently  of  vehicle  typology, 
specifications  and  destination  of  use. 

A  case  study  was  proposed  to  show  the  potential  of  FC  hybrid¬ 
ization  provided  that  suitable  optimal  design  and  energy 
management  are  guaranteed.  The  results  obtained  for  a  bus  shuttle 
are  in  line  with  what  was  previously  addressed  and  tested  by  other 
researchers  active  in  the  field  [38],  thus  confirming  the  potential  of 
the  proposed  methodology. 
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